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^ ■ Abstract 

O 

q^ \ The escape rates of the biaxial single domain spin particles with and without 

^— > ' 

an applied magnetic field are investigated. Using the strict potential field de- 
scription of spin systems developed by Ulyanov and Zaslavskii we obtain new 

q effective Hamiltonians which are considered to be in exact spin-coordinate 

O ' 

correspondence unlike the well studied effective Hamiltonians with the ap- 

/\ ' proximate correspondence. The sharp first-order transition is found in both 



cases. The phase diagram of the transitions depending on the anisotropy con- 
stant and the external field is also given. 
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I. INTRODUCTION 

The decay rate of metastable states or transition rate between degenerate vacua is dom- 
inated at high temperatures by thermal activation, whereas at temperatures close to zero, 
quantum tunneling is relevant. At some critical temperature the transition from the clas- 
sical to the quantum-dominated regime occurs. The transition can be first-order, with a 
discontinuous first derivative of the escape rate, or smooth with only a jump of the second 
derivative in which case it is known as of second-order. 

Based upon the functional-integral approach Affleck JTJ and Larkin and Ovchinnikov |2j 
demonstrated with certain assumptions for the shape of the potential barrier that a second- 
order phase transition from the thermal to the quantum regime takes place at a critical 
temperature T = l//? , where (3$ is the period of small oscillations near the bottom of 
the inverted potential well. Chudnovsky, however, showed that the situation is not generic 
and that the crossover from the thermal to the quantum regime can quite generally be 
the first-order transition || that takes place at T c > To for the case in which the period 
versus energy curve possesses a minimum. Shortly after the observation of Chudnovsky the 
sharp first-order transitions were found theoretically in spin tunneling for two systems. One 
of these is a ferromagnetic bistable large-spin particle j|,[5] described by the Hamiltonian 
H = —DSl — B X S X which is believed to be a good approximation for the molecular magnet 
MnyiAc of spin S = 10, and the other is a biaxial anisotropic model, whose effective mass 
was shown to be position-dependent ||. It was the external field B x (in first model) and 
the anisotropic constant ratio A (in second model) that effect the phase transition of the 
crossover. The same models with the magnetic field applied along alternative axes have 
also been studied, and the corresponding phase diagrams have been given 0,||. A sufficient 
criterion for the first-order transition in the context of tunneling can be obtained by studying 
the Euclidean time period in the neighbourhood of the sphaleron configuration at the peak 
of the potential barrier |9|JT0|]. In the present paper we incorporate the two parameters in a 
single spin tunneling model in order to investigate the dependence of the phase transition 



on these. 

The phenomenon of spin tunneling has attracted considerable attention not only in view 
of the possible experimental test of the tunneling effect for mesoscopic single domain par- 
ticles - in which case it is known as macroscopic quantum tunneling - but also because the 
spin system with an applied field provides various potential shapes and therefore serves as 
a testing ground for theories of instanton induced transitions. The key procedure in deal- 
ing with spin tunneling is to convert the discrete spin system into a continuous one by a 
spin-coordinate correspondence. There are various spin variable techniques which result in 
effective Hamiltonians. It is a long standing question whether the different effective Hamil- 



tonians for a given quantum spin system lead to the same result [|Il],[l^]. Following Ulyanov 
and Zaslavskii we have obtained a new effective Hamiltonian for the spin particle with bi- 
axial anisotropy in addition to the one of a sine-Gordon potential with position dependent 
mass already known in the literature ||[7|,[13].[14 . 



The paper is organized as follows. In Sec. II we give a brief review of the general theory 
of phase transitions of escape rates. Using the effective method of Ref. ]TTJ] we then derive 



in Sec. Ill an alternative effective Hamiltonian for the ferromagnetic particle with a biaxial 
anisotropy without an applied magnetic field. It is shown that the sharp first-order transition 
from the classical to the quantum regime indeed exists in agreement with the observation 
in our previous paper ||. In Sec. IV we then investigate the spin tunneling and phase 
transition with an external field applied along the easy axis which breaks the symmetry and 
makes one of the degenerate vacua metastable. Our conclusions and discussions are given 
in Sec. V. 

II. THE CRITERION FOR THE SHARP FIRST-ORDER PHASE TRANSITION 

OF THE ESCAPE RATE 

At temperature T the escape rate of a particle through a potential barrier can be obtained 
by taking the ensemble average of the tunneling probability T t (E), i. e. 



T(T) = Jr t (E)e-^dE (1) 

where the tunneling probability at a given energy E is defined by 

r t {E) = Ae~ w{E) (2) 

and 

f<l>f( E ) I 

W{E) = 2 / d(f>^2m{(t>)\y{4>) - E] (3) 

J<j>i(E) 

is evaluated from the periodic pseudoparticle (instanton or bounce) trajectories C between 
turning points <fti and <pf. The pseudoparticle trajectory <p c minimizes the Euclidean action 
at the given energy E above the metastable minimum such that 5S((j) c ) = with periodic 
boundary condition 4> c (0) = 4> C (P)- The Euclidean action Se and Lagrangian Ce are 

S E = J drC E = W + (3E, C E = ^m(0)0 2 + V{<j>) (4) 

respectively. Here = difi/dr and t — it denotes Euclidean time. In general the mass 
m(4>) could be position dependent in the context of spin tunneling. The time period {3(E) is 
related to temperature T by (3(E) = ^, as usual. The prefactor A in Eq. (0) results from a 
Gaussian functional integration over small fluctuations around the pseudoparticle trajectory 
(j) c . In the semiclassical approximation the escape rate at temperature T is dominated by 

r(T) ~ e - s ""- (T) , (5) 

where S m i n (T) is the minimum effective Euclidean action which is chosen as the smallest 
value of So and S(T) = Se- Here, So is the thermodynamic action defined by 

S = !3Eo (6) 

with Eq being the barrier height for the pseudoparticle to tunnel through. Generally speak- 
ing, at E = (the bottom of the initial well) the Euler-Lagrange equation leads to the 
vacuum instanton or bounce solution. When < E < Eq (between the bottom and the 
top of the barrier) the trajectory </>(r) shows periodic motion in the barrier region of the 



potential V{4>) which is forbidden for the classical particle. The period of oscillation as a 
function of energy E is given by the following integral 

r<P f (E) \/2m((j)) 

(3(E) = / dtp / WJ . (7) 

There are two independent criteria for the existence of first-order transitions between the 
classical and quantum regimes. The non-monotonic behavior of the oscillation period as a 
function of energy, i.e. the existence of a minimum in the f3 ~ E curve, was once proposed 
as a condition || for first-order phase transitions in quantum mechanical tunneling. The 
corresponding dependence of the action on temperature leads to an abrupt change at some 
critical temperature, at which the first derivative of S m i n (j3) is discontinuous, indicating 
that the crossover from the thermal to the quantum regime is the first-order transition in 
temperature. 



More generally, for a massive particle with position coordinate q, it has been shown flO 
that the existence of a first-order transition leads to the condition 



V"\q a )( gi + f ) + \v""{q s ) + M'(q s )oJ 2 g 2 + M'(q s )( 9l + |) + \M"{q s )u 2 

where 

u 2 M'{q s ) + V'"{q s ) 



<0 (8) 



9i(v) 



±V"{q s 



2M'(q a ) + V"\qs) 
92{UJ) ~ -^M{q s )ul + V"{q s )Y (9) 



and 



, _ V"(q.) 



•* = -* = mj (10) 

and M is the position dependent mass. The subscript s stands for the coordinate of the 
bottom of the well of the inverted potential, i.e., the coordinate of the sphaleron. This 
criterion has been applied to various models studied earlier and the results coincide with 
previous ones. 



III. EFFECTIVE HAMILTONIAN OF THE FERROMAGNETIC PARTICLE 

The model we consider here is that of a nanospin particle which is assumed to have 
a biaxial anisotropy with XOY easy plane and the easy X-axis in the XY-plane, and is 
described by the Hamiltonian 

H = K X S 2 Z + K 2 S 2 , K x > K 2 > (11) 

which has been extensively studied in the context of tunneling from various aspects such as 
ground state tunneling 1 13 JIB], tunneling at finite energy, namely, with the periodic instanton 



14| , and topological quenching of tunneling ||16| , |T7| . Most recently it was shown that this 



model possesses a first-order phase transition from the thermal to the quantum regime ||. 
In all these investigations the quantum spin system of Eq. (|TTD is converted into a potential 
problem by using the conventional spin coherent state technique with approximate spin- 
coordinate correspondence (see Appendix). The effective potential is of the sine-Gordon 
type. 

In the present investigation we reexamine the quantum spin system in terms of a new 
method developed by Ulyanov and Zaslavskii ||11| . The spin operator representation in 
differential form on the basis of spin coherent states is given by the relations 

S z = — i— , S x = S cos ip — sirup— —, S v = Ssinip + cos ip— (12) 

dip dip dip 

The eigenvalue equation 

H\$)=E\$) (13) 

then becomes a second-order differential equation, i.e. 

{Kx - K 2 sin 2 (f)^ + K 2 (s- i) sin 20^ + (E - K 2 S 2 cos 2 - K 2 S sin 2 0) $ = 
where we have shifted the azimuthal angle by | for convenience, cp — <p + \. Following Ref. 



Ill, we use the transformation 



# = $0)0^1 - K 2 sin 2 0)-2 



12 _ x _ ^ 



x= , = F(4>,k), k = A = — - , sin0 = snx (14) 

70 i/l - A sin 2 </>' ^i 



The eigenvalue equation is then transformed to the following effective potential form 

f/ v& rn t 

_ ft _ + ^ (5 + 1) _ $ = M (15) 

where smc, cnx and dnx denote Jacobian elliptic functions. The effective Hamiltonian is 
seen to be 

„2 i 

H = — + U(x), m = , U(x) = K 2 S(S + l)cd 2 x. (16) 

2m 2Ai 

with cdrr = cnx/dnx. We remark here that this derivation, unlike that in previous inves- 
tigations, is exact and without a large s limiting procedure. We also emphasize that this 
Jacobian elliptic potential is of interest on its own and has not been investigated before in 
the context of instanton considerations. The periodic instanton solution leads to an integral 
with finite energy and is obtained as 

x p = sn _1 [/csn(cjr), k] (17) 

The trajectory of one instanton as half of the periodic bounce is shown in Fig. 1 with added 
instanton— anti-instanton pair. We choose s = vTOOO and K\ = l(thus K 2 = A) in the 
diagrams. The period of this periodic configuration is seen to be 

PW = 4L = L lK I P Mk) (18) 

uj{E) yKi \JK 2 s l — EX 

where JC(k) is the complete elliptic integral of the first kind, and 



*"Vn 2 -A' E ' 



to = tu 0] Jl - A w 2 = 4 Kl K 2 s 2 . (19) 

The action calculated along the above trajectory is S = W + (3E with 

W = ^-\K{k)-(l-Xk 2 )U{Xk 2 ,k)} (20) 

XK\ L - 1 
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Here we are particularly interested in the phase transition. To this end we show in Fig. 2 
the shape of the potential for various values of A. The peak of the barrier becomes flatter 
and flatter as A increases. The curve of (3 versus E is given in Fig. 3 and demonstrates the 
obvious first-order phase transition for A > 1/2. Fig. 4 shows the action as a function of 
temperature. 

Next we apply the criterion for the first-order phase transition to the model above. The 
sphaleron is located at x s = 0. Computing the corresponding quantities at the sphaleron 
position, i.e. 

V[x s ] = K 2 s\ V'[x s \ = 0, V"[x s \ = -2K 2 s 2 (l - A), 

V'"[x s ] = 0, V""[x s ] = 8K 2 s 2 (l - A)(l - 2A) (21) 



Eq. (H) becomes 



\v""{x s ) = K 2 s 2 (l - A)(l - 2A) < (22) 



and we regain the critical value of A c = \ at which the first-order transition sets in. We see 
the new effective Hamiltonian with exact spin-coordinate correspondence leads to the same 
results as those of ref. ||. 

However, the physical interpretation for the sharp first-order phase transition is now 
different. In the present case the effective mass is constant and the sharp transition from 
quantum to classical behavior results from a flattening of the peak of the barrier. In Ref. 
H the first-order transition resulted from the position dependence of the mass which makes 
the latter heavier at the top of the barrier. 

IV. THE PHASE TRANSITION WITH AN APPLIED FIELD ALONG THE 

EASY AXIS 

The Hamiltonian with an applied magnetic field h along the easy X-axis is given by 

H = K^l + K 2 S 2 V - gfi B hS z , (23) 



where hb is the Bohr magneton, and g is the spin g-factor which is taken to be 2 here. The 
anisotropy energy associated with this Hamiltonian has two minima: the one on the +X- 
axis which is a metastable state and the other on the —X-axis. Between these two energy 
minima there exists an energy barrier, and the spin escapes from the metastable state either 
by crossing over or by tunneling through the barrier. 

Following the same procedure as in the previous section, the Hamiltonian of Eq. (|23D 
can be mapped onto a point particle problem with effective Hamiltonian 

* = £ + "<*> < 24 > 

where 

m=— — , V(x) = K 2 s 2 (l + — — )sn 2 (x, k) + K 2 s 2 a (cn(x, k)dn(x, k) - l) (25) 

where the metastable minimum of the potential has been shifted to x = for convenience. 
Now the effective mass is a constant. The potential is more like that of an inverted double- 
well potential with the sphaleron position x s = cd _1 (^, k) and barrier height 



« 2 A, l-f 2 (l + f)(f~ 1) 

— )ttaz + K2S a — fr^ 



E = K 2 s\l + —)—^t + K 2 s*a K 2 >\ 2 (26) 



i 

Under the barrier a bounce configuration exists as shown in Fig. 5. We redraw the period 
diagram as a function of energy for the same parameter (A = 0.9, a = 1) in Fig. 6 and 
obtain the first-order transition from thermal activation to quantum tunneling as shown in 
Fig. 7. 

Applying the phase transition criterion to this model, we obtain 

a 2 
V'(x s ) = 0, V"(x s ) = -2K 2 s 2 (l - -)(1 - A), 



V"\x s ) = -3K 2 s 2 a(l - A) 



1 oP_ 

1 4 



\i-4» 



i " 



= 2K 2 s 2 (l A) / 2 _ 2A)ft4 _ (7a2 _ 22A + 7)ft2 _ 32A + x 
4 — a 2 \ v ' 

and the frequency u 2 is 



2 2 



V"{x s ) 



J — ^s ~ 

m 



-2K x V"{x B ) (28) 



The expressions for g\ and gi are now found to be 

|«(1-A) 



9\ 



l-f)(l-SjA) 



92 = ; 8 29 

^(i-x)! 1 -^) 

The critical line of the two parameters for the first order transition requires 



2(2 + o^ 

+ a 2 



A>^^ (30) 



The corresponding phase diagram is shown in Fig. 8. From the diagram we observe several 
interesting features. First, the classical-quantum phase transition shows both the first-order 
(region I) and the second-order (region II) transition domains. We see that there is only a 
second-order transition for A < 0.5. For materials with A larger than 0.5 we can see that 
the order of the phase transition changes from first to second as a increases and the phase 
boundary changes with A up to 1. 



An alternative effective Hamiltonian [18,19] with the conventional application of the spin 



coherent state technique as that for the biaxial anisotropy spin particle without the applied 
magnetic field || has also been investigated. However, the effective Hamiltonian with ap- 
proximate spin-coordinate correspondence gives rise to a result ]19] for the phase transition 
which differs from that of the Hamiltonian of Eq. (21) with exact spin-coordinate corre- 



spondence. The center of the position dependent mass in the Hamiltonian with approximate 
spin-coordinate correspondence in Refs. |T8| , P^| does not coincide with the position of the 
pseudoparticle, namely, the bounce here. We, therefore, believe that the instanton concept 
fails in this case. A detailed analysis will be reported elsewhere. 

V. CONCLUSION 

In summary, we investigated the biaxial spin model with applied magnetic field along 
the easy X-axis with the effective potential methods developed in Ref. [|ll |. The resulting 



10 



periodic instanton solutions with potential in terms of Jacobian elliptic functions and the 
corresponding phase transitions have been studied for the first time. For the quantum spin 
particles with anisotropic constant A > | and above the critical value of magnetic field a, 
the energy dependence of the oscillation period shows the obvious first-order transition from 
thermal activation to thermally assisted quantum tunneling. Applying the criterion for the 
first-order transition, we obtained the phase diagram which exhibits two domains separated 
by a critical line, indicating the first-order and second-order transitions respectively. 

Appendix 

In the conventional application of the spin coherent state technique, two canonical vari- 
ables, and p = s cos 9 are adopted with the usual quantization 

[<P,p]=t (Al) 

We show in the following that the spin-coordinate correspondence is only approximate up 
to order 0(s~ 3 ). 

From the relation between the spin operators and the polar coordinate angles 

S x = s sin 9 cos 0, S y = s sin 9 sin 0, S z = s cos 9 (A2) 

the usual commutation relation of spin operators reads 

[S x , Sy] = s [sin 9 cos 0, sin 9 sin 0] = s sin 9 [cos 0, sin 9] sin + s sin 9 [sin 9, sin 0] cos 

(A3) 
Using Eq. (|A1|) , one can prove the following relation 

[sin 9, cos 0] = A + cos + iA_ sin 0, [sin 9, sin 0] = A + sin — iA^ cos (A4) 

with 

A + = - (sfl - (cos6> + a) 2 + sjl - (cos 9 - a) 2 
1 



A_ = - ( a/1 - (cosfl + a) 2 - A /l - (costf - c^ 2 
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where a = 1/s. Substituting ( |A4|) into Eq. (|A3|) , one has 

[S X ,S V ] = s 2 (-i) sinO A_ = icos6s 2 a + 0(a 3 ) (A5) 



i.e. 



[S x ,S y ]=iS s + 0{s- 3 ) (A6) 

which implies that the usual commutation relation holds only in the large spin limit. 
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Figure Captions: 

Fig. 1: The effective potential and the corresponding periodic instanton configurations. 

Fig. 2: The potentials for different values of A. 

Fig. 3: The oscillation period as a function of energy at A = 0.9. 

Fig. 4: The action as a function of temperature indicating the first-order transition at 

A = 0.9. 

Fig. 5: The effective potential in the case with external magnetic field. 

Fig. 6: The oscillation period as a function of energy at A = 0.9 and a = 1. 

Fig. 7: The action as a function of temperature showing the first-order transition in the 

case with external field. 

Fig. 8: The phase diagram for the orders of phase transitions in the (A, a) plane. Region I: 

the first-order domain, region II: the second-order domain. 
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